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We present a framework to compute the responses of hadron masses to the chemical potential in 
lattice QCD simulations. As a first trial, the screening mass of the pseudoscalar meson and its first 
and second responses are evaluated. We present results on a 16 x 8^ x 4 lattice with two flavors of 
staggered quarks below and above Tc- The responses to both the isoscalar and isovector chemical 
potentials are obtained. They show different behavior in the low and the high temperature phases, 
which may be explained as a consequence of chiral symmetry breaking and restoration, respectively. 
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I. INTRODUCTION 

It is vifell known that studying finite density QCD 
through lattice simulations is a very hard problem. The 
fermionic determinant at finite chemical potential is com- 
plex, and gives an oscillating behavior in quantum aver- 
ages which makes simulations very inefficient. Since the 
naive quenched approximation at finite chemical poten- 
tial leads to an essentially different world m, the use of 
dynamical fermions would be essential to extract the rel- 
evant physics. In spite of these difficulties, the study of 
hadrons in a finite baryonic environment is quite impor- 
tant P, 0, El , in view of recent experimental developments 
and of the theoretical interest in the phase structure of 
QCD. Search for the quark gluon plasma phase in high 
energy heavy ion collision experiments requires theoret- 
ical understanding of hadronic properties at finite tem- 
perature and density H. Moreover, some experimental 
results can be interpreted by assuming a shift in the mass 
and the width of the p meson, induced by the dense nu- 
clear medium even below the deconfinement transition 



such as masses, in the vicinity of zero chemical potential 
at finite temperature, and explore their changes through 
the response to the chemical potential at /i = 0. This al- 
lows to perform the numerical simulations with standard 
methods. The Taylor expansion in /i has also been used, 
and its properties discussed, in ||]. There, the behavior 
of observables, measured by standard methods as a func- 
tion of an imaginary chemical potential /i, is fitted by a 
Taylor series amenable to analytic continuation to real 
/i (assuming a large enough analyticity domain). Here 
we measure the Taylor coefficients directly, in a single 
simulation, at /z = 0, by calculating the derivatives of 
the relevant observables. Although the Taylor expansion 
cannot reproduce the non-analyticity inherent to a phase 
transition, it may suffice for observing the rounded, ana- 
lytic behavior indicative of a phase transition in a finite 
volume, and for applying finite size scaling to probe the 
transition in the infinite volume limit. Our preliminary 
results have been reported in |11| . 



In this paper we propose a new technique to investi- 
gate non-zero chemical potential using lattice QCD sim- 
ulations. There are several approaches to circumvent the 
difficulty of studying a finite chemical potential system, 
and they seem successful to a limited extent [0, ||, ||. 
In particular, the study of baryon number susceptibility 
at zero baryon density has reported an abrupt jump at 
the transition temperature pl]| . There is in fact much 
interesting physical information which can be extracted 
from the behavior of a system at small chemical poten- 
tial. Our strategy is to expand the hadronic quantities, 



This paper is organized as follows. In the next section, 
we develop the basic formulae to evaluate the first and the 
second responses of hadron screening masses with respect 
to both of the isoscalar and the isovector chemical poten- 
tials. In Section [I], we report on Monte Carlo simula- 
tions performed with two dynamical flavors of staggered 
quarks. We obtain the responses of the pseudoscalar me- 
son mass, and discuss some implications of our results. 



Section IV presents our conclusions. 



II. CHEMICAL POTENTIAL RESPONSE OF 
HADRON MASSES AT FINITE TEMPERATURE 

This section develops the basic framework to observe 
the response of hadron masses with respect to the chem- 
ical potential. At fixed temperature T and bare quark 
masses, the screening mass of a hadron is expanded in 
the form: 
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Assume that the spatial hadron correlator C{x) is domi- 
nated by a single pole contribution, 



Cix) = J2{H{x,y,zA)HiO,0,0,0)^) 



(2) 



where M = aM and a; = x/a. L^ is the lattice size in 
the x-direction. Generalization to a multi-pole situation 
is straightforward. In the case of a local hadron operator. 



A is represented as ^ = j/2M with the residue of the 
propagator, 7. We employ a smeared hadron operator 
H{x, y, z, t) to enhance the overlap with the ground state, 
and hence A depends on the choice of sources. With a 
fixed type of smearing function, however, the behavior of 
A as a function of the chemical potential provides infor- 
mation on the coupling to the medium. 

We take the first and second derivatives of the hadron 
correlator with respect to jl = afi — ii/{NiT), where /i 
is the chemical potential. While we here write down ex- 
pressions which depend on a single chemical potential, it 
is straightforward to generalize to the case involving two 
flavors corresponding to u and d quarks. Later we intro- 
duce the isovector and the isoscalar chemical potential 
by setting /z„ and /id appropriately. The first and the 
second derivatives of C{x) with respect to the chemical 
potential lead to 



1 c 


lC{x) 


1 dA 






C(x) 


dfl 


A dfl 




dM 
dfi 


{(- 


LA , 
tanh 

2 J 


M\x- 

. V 


-t)] 



T 



(3) 



and 



J 



1 d^Cjx) 
C{x) djj? 



1 d'^A 



tanh 




M\x- 



L, 



tanh 



M X 




(4) 



C(x) and the derivatives of C{x) are calculated from 
lattice simulations. Then, using the right-hand side of 
Eqs. (^) and (^), the first and the second responses of 
hadron masses and couplings are determined. 

The next step is to relate the derivatives of the cor- 
relator to quantities which are measured in the lattice 
simulations. In this work, we consider the flavor non- 
singlet mesons in QCD with two flavors. The hadron 
correlator is then given by 



where G is the meson propagator 

G - Tr [PifiuUTP{pd)onT^] 



(5) 



(6) 



Here P{fi) is the quark propagator at finite chemical po- 
tential, and r is the Dirac matrix which specifies the 
spin of the meson. The quark propagator is related to 
the Dirac operator D[U; p] in the background gauge field 
configuration U as 



P{ii)^D[U;(ir\ 



(7) 



The expectation value of operator O, (O), stands for 
J[dU]OAe-^<^ 
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where Sq is the gluonic action and A is the fermion de- 
terminant, 

A = detD[U;fLu]detD[U;fid]. (9) 

Then, the first and the second derivatives are represented 
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where the dotted O and O denote the first and the second 
derivatives of the operator O with respect to fi. 

At zero chemical potential, we have simpler expres- 
sions, since 
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Eq. (|2|) corresponds to the fact that the average baryon 
number density vanishes at /i = 0. Actually, we see that 
d{(iet D) / dfi = Tr[DD~^] detZ? is anti-hermitian at jj, — 
0: 
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This means that d{detD)/d[l changes sign under the 
transformation U ^ U^ . Since the measure and the glu- 
onic action are invariant under this transformation, its 
expectation value vanish es p^ ]. Thus, at zero chemical 
potential, Eqs. ( p^ and (|ll|)~turn into 
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We investigate derivatives with respect to both the 
isoscalar and the isovector types of chemical potential. 
The isoscalar chemical potential is conjugate to the total 
quark density. In this paper we study the response to the 
isoscalar chemical potential by setting 



Ms — A^« — Md, 
and for the isovector case 

Av = Au = ^Ad- 



(16) 
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Note that Son and Stephanov proposed a model corre- 
sponding to the isovector case as a good test bed for 
chemical potential effects in QCD 13 . An advantage of 
setting Au = —fj-d is that the fermion determinant is pos- 
itive, and therefore the problem becomes tractable with 
standard lattice techniques. Here we do not make use of 
this advantage and still study the dependence on Ay by 
performing a Taylor expansion around Av = 0- Future 



simulations at non-zero Ay would also constitute a good 
test of the performance of our approach. 

Our simulations are performed with Nf = 2 dynam- 
ical quarks of the staggered fermion type. The fermion 
operator and its derivatives are 

D[U;p]n,m = maSn^rn 
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where a and t are unit vectors pointing along spatial and 
temporal directions, respectively. 

Taking into account the four-fold degeneracy of the 
staggered fermion operator, the determinant factor A of 
Nf = 2 fermions is then 

A = exp QTrlni?[C/;Au] + ^Tt \nD[U ; fid]) . (21) 

Explicit formulae for the isoscalar and the isovector cases 
are presented in Appendices A and B respectively. Ap- 
pendix C gives the specific forms for staggered quarks. 



III. NUMERICAL RESULTS 
A. Lattice setup and meson screening masses 

In this section, we describe the results of Monte Carlo 
(MC) simulations for the response of meson screening 
masses with respect to the chemical potentials below and 
above the chiral phase transition. Numerical simulations 
are performed on lattices of size 16 x 8^ x 4 with the 
standard Wilson gauge action and with two dynamical 
flavors of staggered quarks. The configurations are gen- 
erated with the R-algorithm with molecular dynamics 
step 5 = 0.01 and trajectory length of 50 steps. For 
Nt — 4, the critical coupling (3c of the finite tempera- 
ture phase transition has been determined as /3c = 5.271 
for ma = 0.0125 and /3c = 5.288 for ma = 0.025, re- 
spectively [n4l n3] . We perform numerical simulations at 
ma = 0.0125, 0.017 and 0.025. The values of /3 are listed 
in Table Q together with the number of configurations and 
rough estimates of the corresponding temperatures. To 
evaluate the traces of the various fermionic operators, the 
Z2 noise method [|6[ is used, with 200 noise vectors. The 
real part of the average Polyakov loop, as shown in Fig- 
ure 0, exhibits a rapid increase between about (3 = 5.27 
and 5.28, as expected. 



TABLE I: Parameters of the simulations. Nconf denotes the 
number of configurations on which the correlators are mea- 
sured. 
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FIG. 1: Average value of the real part of Polyakov loop as a 
function of /3. 



Meson correlators are measured on configurations sep- 
arated by 20 trajectories. We use a corner-type wall 
source ||lj after Coulomb gauge fixing in each (y, z, t)- 
hjperplane. We determine the meson screening masses 
M by fitting the MC data to the single exponential 
form, Eq. (0), with the fitting range x = 1-15. As 
shown in Figure g, the fit is successfully applied for the 
pseudoscalar(PS) meson channel. The extracted PS me- 
son screening mass is listed in Table ||. Below Tc, the 
present lattice extent in the shorter spatial directions is 
not sufficiently large to extract reliably the PS meson 
mass. Since it is known that the finite size effects on 
hadron masses are rather severe with dynamical stag- 
gered fermions [|l8|, |l^, one should use lattices with L-m-^ 
larger than ~ 5 to avoid finite size effects, and hence 
larger lattice sizes are obviously required for a quanti- 
tative analysis. Goals of the present study are, how- 
ever, to show the applicability of the method to extract 
the responses of meson mass to the chemical potentials, 
and to observe the qualitatively different behavior below 
and above Tc giving us hints for future detailed studies. 



TABLE II: Responses of the PS meson to the isoscalar chem- 
ical potential fis. 
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1 d^A 
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0.0125 5.26 0.2956( 2) 
5.34 0.7513(11) 



-1.4(20) 0.17(35) -0.8(23) 
-4.23(49) 5.39(10) 2.95(51) 



0.0170 5.26 0.3506( 2) 
5.34 0.7421(35) 



-1.5(14) 0.30(26) -0.6(16) 
-3.68(75) 5.82(16) 4.16(78) 



0.0250 5.20 0.4061( 2) -0.4(14) 0.16(26) 0.0(16) 

5.26 0.4218( 2) -0.8(11) 0.29(20) -0.1(12) 

5.32 0.6926(11) -4.65(91) 5.17(20) 2.82(96) 

5.34 0.7534( 7) -3.17(41) 4.43( 8) 2.71(42) 
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FIG. 2: Pseudoscalar correlator and several contributions to 
its responses to fis,v (see Eq. (E2)) at /3 = 5.34 and ma = 
0.025. The correlator is fitted to the single pole form, Eq. (bl). 



Therefore we present also the results below Tc in spite 
of the large systematic uncertainty. The screening mass 
of the PS meson at ma = 0.025 as a function of T/Tc is 
shown in Figure 0. 

In the case of the vector meson correlator, although 
the procedure is in principle also applicable, the statis- 
tical fluctuations prevent us from reaching any relevant 
conclusion at this stage, and therefore we do not report 
results for the vector channel in the following. 



B. Responses of the PS meson mass to the 
isoscalar chemical potential 

Several terms contribute to the derivatives of the me- 
son correlator with respect to the chemical potential: see 
Eqs. (B5), (|C3|), and (p4[) in Appendices. As represen- 
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FIG. 3; Screening mass of the PS meson in lattice units 
at ma = 0.025. f3c is taken from urn. Values of T/Tc are 
estimated by using the two-loop /3-function. 



tative terms, in Figure ^ we show 
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(22) 



All of them are determined with reasonable accuracy. 

Let us start with the first response of the PS correlator 
to the chemical potential. Note that the first derivative 
with respect to the isosc alar chemical potential is identi- 
cally zero (see Eq. ( AlC ) in Appendix A) . For the isovec- 
tor chemical potential, Figure H shows C^^dC/dfiy at 
/3 = 5.26 and 5.34. At both temperatures, the values are 
very small. This is consistent with a comparative study in 
the Nambu-Jona-Lasinio model, which also predicts very 
small responses around the critical temperature JE1|. 

We then consider the second responses. Figure |^hows 
C-^dPC/dfi'^ at /? = 5.26 (below T^) and 5.34 (above TJ 
for the isoscalar and isovector chemical potentials. The 
solid curves represent the fits to Eq. (0), after fitting C{x) 
to Eq. d). 

We determine the responses of the meson mass with re- 
spect to the chemical potentials as follows. As described 
in the last subsection, the meson correlator is first fitted 
to Eq. (g) to extract the screening mass, M. Then, we 
fit the MC results for the derivatives of the meson cor- 
relator to Eqs. (H) and (||), substituting the determined 
value of M. The derivatives of mass and coupling are 
then obtained as fitting parameters. Note that for fis we 
omit the fitting step to Eq. (||), since the first response 
is strictly zero. 
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FIG. 4: The first response of the PS meson correlator at 
P = 5.26 (T < Tc, top figure) and (i = 5.34 (T > Tc, bottom). 
The quark mass is ma — 0.025. 



Results of the PS meson response to the isoscalar chem- 
ical potential are summarized in Table O. In the low 
temperature phase, the dependence of the mass on fis is 
small. This behavior is to be expected, since, below the 
critical temperature and in the vicinity of zero fis, the 
PS meson is still a Goldstone boson. In fact, the extrap- 
olated value to the chiral limit of the isoscalar response 
is consistent with zero, as shown in Figure pi This is in 
contrast with the behavior above Tc, where d^M'^ /d[l^ 
seems to remain finite even in the chiral limit. In addi- 
tion, our results suggest that the response of the coupling 
is small below T^- 

Above Tc, we first note that the correlator and its re- 
sponse are still well fitted by the single pole formulae, 
Eqs. (Ij-Q). The screening masses are manifestly larger 
than those below Tc. This is consistent with previous 
work [|0|. As pointed out above, the response of the 
mass above Tc becomes large, as a reflection of the fact 
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FIG. 5: The second responses of the PS meson correlator 
at 13 ^ 5.26 and ma = 0.025 (T < T^, top), and at /3 = 5.34 
and ma = 0.0125 (T > Tc, bottom). The curves are fits to 
Eq. ffl). 



that the pion is no longer a Goldstone boson and as an 
indication of chiral symmetry restoration. We also find 
that the response of the coupling 7 increases above Tc, as 
shown in Figure & which may indicate a larger overlap of 
the wall-source operator with free quarks in the plasma 
phase. Note that {(P^/djl^)/^ does not depend on the 
choice of the source normalization. 



Responses of the PS meson mass to the 
isovector chemical potential 

Results for the isovector chemical potential are sum- 
In the presence of the isovector 
and 7r~ may have different masses. 



FIG. 



marized in Table III 



M^ (top) and (fM^/djX^ (bottom) for the PS meson 
versus ma. f3 is 5.26 (triangles) and 5.34 (circles). Extrapo- 
lation to ma = is also shown. 



chemical potential, tt 

Here we consider the n'^ ( ud ) meson as shown in Eq 



In contrast to the case of the isoscalar chemical poten- 
tial, the second order response of the mass is significantly 
large in the low temperature phase, and decreases in mag- 
nitude above Tc. The difference between the isoscalar and 
isovector chemical potentials is illustrated in Figure 0. 
The response of 7 also shows different behaviors in the 
confined and deconfined phases, illustrated in Figure ^. 

These features are manifest even for a small quark mass 
parameter. Note that the isovector potential explicitly 
breaks the u-d symmetry, even if the two quarks have 
equal masses. The phase structure in the (T, \fiv\) plane 
has been studied by Son and Stephanov fls] . The original 
SU{2)L+ji symmetry at non-zero quark mass and zero 
chemical potential is broken down to U{1)l+r. At zero 
T and for \^v\ larger than the mass of the pseudoscalar, 
the system is in a different phase than at ^v = 0- The 
ground state is a pion condensate and there is one mass- 
less Goldstone boson associated with the spontaneous 
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FIG. 7: Second responses (fM/dfis and (fM/djly of the PS 
meson mass at ma — 0.025. T/Tc is estimated as in Figure ^ 



breaking of the U{1)l+r symmetry. For \^v\ = mps, 
the critical temperatm'e is T = 0. At sufRciently high 
temperatm'e, the condensate melts and the symmetry is 
restored. Due to the presence of the phase boundary, we 
do not expect to be able to reach the condensed phase 
by Taylor expanding around fj,v = 0. We can, however, 
hope to get some hints about the presence of the phase 
boundary while keeping [/xy| < mps- In this case, the 
system is in the same ground state as for zero chemical 
potential, there are no exact Goldstone modes and the 
three pions are massive. 

An interesting point in this respect is that the second 
derivative of the mass is negative in the low temperature 
phase, in marked contrast with what happened for the 
isoscalar potential. The mass tends to decrease under the 
influence of the isovector chemical potential, reflecting 
the fact that for low temperature and chemical potential 
above the pion mass, a Goldstone mode appears |^, [l3[ . 
This is more clearly shown by an expansion as in Eq. ^). 
At (3 = 5.26 and ma — 0.017, the data suggest 
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(23) 



The coefficient of the linear term is consistent with zero. 
Notice also in Table III that the lighter the quark mass, 
the stronger the response, a possible indication that for 
lighter pions the phase boundary is closer to the zero 
chemical potential axis, as suggested in Q. 

In the high temperature phase, the dependence of 
the masses on ^y decreases. Since the PS meson be- 
comes heavier, the phase boundary to the pion conden- 
sate phase is farther away from the /iy — axis. The 




FIG. 8: Second responses of the couplings, 'y~^d?^/dji%, 
j'-^d^'y/djly, for the PS meson at ma = 0.025. T/Tc is esti- 
mated as in Figure H. 



weaker responses may be understood from this point of 
view. 



IV. CONCLUSIONS 

In this work, we have developed a framework to study 
the response of hadrons to the chemical potential. It 
is based on Taylor expanding hadronic quantities about 
fi = 0. We show the first results of the first and second 
derivatives of the PS meson screening mass with respect 
to fj,. As shown in the previous section, the second or- 
der responses are sizable and reveal several characteristic 
features. The behavior of the responses for the PS me- 
son seems to have close contact to the chiral symmetry 
restoration. For the isoscalar chemical potential iis, the 
dependence of the PS mass on /ig in the chiral limit is 
consistent with zero, reflecting the fact that at low tem- 
perature and small /ig the pion is still a Goldstone boson. 
For the isovector chemical potential, we show features 
consistent with the phase structure proposed by Son and 
Stephanov [|l3|. The wJ PS meson mass tends to decrease 
as a function of fiy at a much stronger rate in the low 
temperature phase. 

It is notable that a single hadron pole gives a good 
description for the response as well as for the correlator 
even at /3 = 5.34 {T/Tc ~ 1-1) in the PS channel. 

Since the present study is a first trial, our simulations 
have been performed on a rather small lattice. Since the 
difference between the dynamics of A^t = 4 and Nt = 6 
lattices is sizable |2Q] and since the spatial volume should 
be large enough to accommodate the hadrons, further 
investigations on larger lattices are indispensable. The 
chemical potential responses of other quantities, such as 
the nucleon mass and the quark condensate, are impor- 
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TABLE III: Responses of the PS meson correlator to the isovector chemical potential p.v- 



ma 


/3 


1 dA 
A dfi 


dM 
dfi 


1 d'^A 
A dp.'^ 


d^M 
dp,-' 


1 d^^ 
7dA^ 


0.0125 


5.26 
5.34 


0.0029(57) 
0.0047(93) 


-0.0001(12) 
0.0006(21) 


47.46(71) 
2.32(64) 


-12.93(43) 
-1.32(32) 


3.7(16) 
0.56(77) 


0.0170 


5.26 
5.34 


-0.0081(48) 
0.0000(82) 


-0.0005(10) 
-0.0012(19) 


33.52(61) 
2.74(64) 


-10.49(33) 
-1.48(32) 


3.6(11) 

0.75(78) 



5.20 


0.0062(39) 


0.0016( 8) 


25.24(46) 


-9.10(23) 


2.84(74) 


5.26 


-0.0080(37) 


0.0007( 8) 


23.22(46) 


-8.64(23) 


2.72(71) 


5.32 


-0.0054(64) 


-0.0020(14) 


4.04(75) 


-2.14(38) 


0.95(93) 


5.34 


0.000(6) 


0.000(1) 


2.99(53) 


-1.51(26) 


0.99(60) 



tant issues. Exploratory studies on these subjects are in is used. 



progress. 



The first derivative at /is = is 
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APPENDIX A: FORMULAE FOR THE 

ISOSCALAR CHEMICAL POTENTIAL 

RESPONSE 

The quark propagator P satisfies the relation 

P{Pd)on^lbP\-ft-d)nOlb, (Al) 

SO that the meson correlator G is given by 



G = Tr 



P(A5)nor75P(-As)io75r^ 



(A2) 



where "Tr" means the trace over spinor and color indices. 
Each propagator is expanded as 



^(A) 



P - fiPDP + —{2PDPDP - PiJP) 



-O(A'), 



(A3) 



A' 



Pi-fi) = P + fiPDP+^{2PDPDP - PDP) 



-O(A^), 



(A4) 



where P and D are the quark propagator and the Dirac 
operator at zero chemical potential, respectively, and the 
relation 



and the second derivative at fis = is obtained as 

G = 4ReTr [{PDPDPUT-f^Plol^r^ 
-2ReTr [{PDPUrjzPlo^Br^ 
-2Tr [{PDPUrToiPDP)iol5T^\ ■ (A7) 

Let us calculate the derivatives of the fermionic deter- 
minant A. Using the following equations, 



dfl 
d^ 



detD = Tr[DP]detD, 



—r^detD = (tt[DP] - TtIDPDP] + Tr[bpf\ 



dfi 



we have 



X detD, 



(AS) 



2Tr 



— = 2Tr 

A 



DP 
DP 



-2Tr 



DPDP 



4Tr 



DP 



.(A9) 



Combining Eqs. (gj), l\A§, (jATl), l\Ag) and (|Ag|), we 
have 



dp, 



Re(iJ(n)iJ(0)t) = 0, 



(AlO) 



P = -PDP 



(A5) 



and 



dfi'^ 



Re{H{n)H{Oy) = 4Re(Tr {PDPDP)norj5PLl5r 



2Re(Tr 



{PDPUrj^Ploi^r^ 



-2Re (Tt 



{PDPUrj5{PDP)l,j,r^ 



8(ImTr 



{PDPUr-f^PLi^r^ 



IniTr 



DP 



+2Re| ^Tr [p„or75^io75rt] (Tr[DP] - Ti[DPDP] + 2Tt[DP] 



Tr 



PnoTlsPlol^r^] ) (Tr[Z)P] - Tt[DPDP] + 2Tr[DPf) } 



(All) 



APPENDIX B: FORMULAE FOR THE 

ISOVECTOR CHEMICAL POTENTIAL 

RESPONSE 

Next, we consider responses to the isovector chemical 
potential, Eq. (17). In this case, the first derivative of A 
vanishes, 



dfiv 



{det D[U;ilv] det D[U ; - fiy]) 



= 0, 



/iv=0 



and the second derivative is obtained as 



— = 2Tr[7:>P] - 2Tr[DPDP]. 
Similarly, derivatives of G are calculated as 



G = -2ReTr 



and 



G 



4ReTr 
-2RcTr 



(pi)P)„or75i^,!o75r^ 



{PDPDPUrToPloTo^^ 

{PDPUr^^Ploi^r^ 



(Bl) 
(B2) 

(B3) 



-2Tr 



{PDPUr^5{PDP)ioi5r^] . (B4) 



Resulting expressions for the first and second responses 
to the isovector chemical potentials are 



—Re(H{n)H{0)^) 



-2ReTr 



(PI)P)„or75Pio75r^l , (B5) 



and 



—Re{H{n)H{0)^) 



4Re(Tr 
-2Re (tt 
+2Rc (tt 



{PDPDPUTj^Plol^r^ 

(pz)p)„or75Pio75r^ 
(pi)P)„or75(Pi>P)Io75r^ 



{C 



-2Rc <^ ( Tr 



PnoT-f5Piol^T^ 



Tr 



PnoT^5Plo^5T^ 



Tt[DP] - Tt[DPDP] 
Tt[DP] -Tt[DPDP]\\ 



(B6) 



APPENDIX C: RESPONSES FOR STAGGERED 
FERMIONS 



For staggered fermions, the determinant factor A is 
given by Eq. (EH), and this leads to 



A 
A 



1 



Tt[DP], 



(CI) 



I = ^Tt[DP] - ^Tr[DPDP] + ^Tt[DP]\ (C2) 
The final expressions for the second responses are 



dfi 



^Re{Hin)H{Oy) = 4Rc(Tr {POPOPUrj^PlolsT^ 



2Re ( Tr 



(PDP)„or75P^o75rt' 



-2Re(Tr 



{PDPUrTo{PDP)lo75r^ 



2(ImTr 



(PI)P)„or75n!o75r^ lmTv[DP] 



-^^° 



Tr 

-(Tr 



PnoT-f^Plol^r^l [Tt[DP] - Tt[DPDP] + -Tr[i>P]2 
Pnori5Plol5r^] ) /Tr[Z)P] - Tr[bPDP] + ^TiiDPf 



(C3) 



10 



for the isoscalar chemical potential, and 



— Re(iJ(n)i/(0)t) = 4Re (Tr [(PI)PI)P)„or75Pio75rt 



2Re (Tt 



(PZ)P)„or75P^o75rt 



+2Re/Tr 



{PDPUrj,{PDP)ioi5r^ 



Tr 

- ( Tr 



Pnor75Pio75rt] (Tt[DP] - Tr[i)PI)P]' 
PnoTj^Plol^T^] ) (Tr[i)P] - Tr[I)PI)P]) | 



(C4) 



for the isovector chemical potential. 
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